Static and Radiating Solutions of Lovelock Gravity in the 
Presence of a Perfect Fluid 



M. H. Dehghani^'^* and N. Farhangkhah^ 
^ Physics Department and Biruni Observatory, 
College of Sciences, Shiraz University, Shiraz 71454, Iran 
^ Research Institute for Astrophysics and Astronomy of Maragha (RIAAM), Maragha, Iran 

Abstract 

We present a general solution of third order Lovelock gravity in the presence of a specific type 
II perfect fluid. This solution for linear equation of state, p = w{p — 43) contains all the known 
solutions of third order Lovelock gravity in the literature and some new static and radiating solu- 
tions for different values of w and B. Specially, we consider the properties of static and radiating 
solutions for w = and w = {n — 2)~^ with B = Q and B ^ 0. These solutions are asymptotically 
flat for B = 0, while they are asymptotically (anti)-de Sitter for B ^ 0. The new static solutions 
for these choices of B and w present black holes with one or two horizons, extreme black holes or 
naked singularities provided the parameters of the solutions are chosen suitable. The static solu- 
tion with w = and vanishing geometrical mass (m = 0) may present a black hole with two inner 
and outer horizons. This is a peculiar feature of the third order Lovelock gravity, which does not 
occur in lower order Lovelock gravity. We also, investigate the properties of radiating solutions for 
these values of B and w, and compare the singularity strengths of them with the known radiating 
solutions of third order Lovelock gravity. 



* email address: mhd@shirazu.ac.ir 



1 



I. INTRODUCTION 



The dynamics of our universe at the classical level may be described by a gravitational 
field equation in which the contribution of energy content of the universe is represented by an 
energy momentum tensor appearing on the right-hand side (RHS) of this equation. The left- 
hand side (LHS) represents pure geometry given by the curvature of spacetime. Gravitational 
equations in their original form with energy-momentum tensor of normal matter cannot 
lead to an accelerated expanding universe, while high-precision observational data have 
confirmed with startling evidence that the universe is undergoing a phase of accelerated 
expansion. There are, then, two ways to obtain accelerated expansion of the universe. 
The first way is by modifying the RHS and supplementing energy-momentum tensor by 
dark energy component. The second way, which one may supersede the dark energy, is 
the generalization of the LHS of the field equation. A natural generalization of LHS of 
the field equation with the assumption of Einstein-that it is the most general symmetric 
conserved tensor containing no more than two derivatives of the metric- is Lovelock gravity 

H. .,.,..„,....,.>.„„..„.._., 

string theory at the classical level [2] • Asymptotically fiat solution of second order Lovelock 



gravity has been introduced in [3|, while this kind of solution for third order Lovelock gravity 
has been introduced in [4]. In recent years, solutions of Lovelock gravity in the presence of 
more general matter distribution have been investigated. Here, we want to introduce the 
exact solutions of third order Lovelock gravity in the presence of a perfect fiuid with linear 
equation of state and investigate their properties. This is motivated by the fact that third 
order Lovelock gravity has some peculiar features which do not occur in lower order Lovelock 
gravity. Recently, one of us has shown that a topological black hole of third order Lovelock 
gravity has an unstable phase which does not occur in the lower order Lovelock gravity j^. 
The black hole solutions of this theory and their thermodynamics have been considered in 
Gj, while radiating solution of this theory in the presence of a null fiuid has been considered 
in 7|. Recently, the generalized Vaidya spacetime in Lovelock gravity has been investigated 
8|]. Also, some efforts have been done in the computation of the conserved quantities of 



m 



third order Lovelock gravity through the use of the counterterm method 

The outline of this work is as follows. In Sec. [Tll we give a brief review of field equations 
and introduce a general solution of third order Lovelock gravity in the presence of a specific 
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type II perfect fluid. Section [TTTl is devoted to the presentation of the static solutions for a 
type I perfect fluid with hnear equation of state. In Sec. HVl we generalize these solutions 
to the case of radiating solutions, and study the properties of them. Finally, we give some 
concluding remarks. 



II. GENERAL SOLUTIONS 



The field equation of third order Lovelock gravity in the absence of cosmological constant 
may be written as ^ 



i=l 



where a^'s are Lovelock coefficients, T^i, is the energy-momentum tensor of matter field, 
£(1) = R, £(2) = R^^^sR^''^^ - ^R^uR^" + R\ and 

£(3) ^ 2R^^''-^R^^p,R'^^^ + SR'^'^^^R'^^^R^^, + 24R^'"^^R^^,pRP^ 

+mR^'"'^R,,^, + 2AR>''"'''R,^R,, + IQR>'''R,,R\ - URR^^'R^, + R' (2) 

are first, second and third order Lovelock Lagrangian, respectively. In Eq. ([1]), H^f}u is just 
the Ricci tensor, and 7i|S and H^^^^^ are given as 

= '^{RfianrRJ"^^ — '^R/ipuaR^'' " "^R/iaR^u + RR/J-u), (3) 



J^a-KTpJ^up "T on ufjpn, T^n, ^ -\- on, y^-^^n ^^n ^ 

-\-ARj^ RfjI^^pR^^j. ^Ry RcTKTpR^ n 4:R ^ RcrKTflRup ^RR^ ^RrpK/l 

+AR'^^RpT-Ry^ — SK^ ^R^-pRP p + ARR^pR''^ — R^R^p). (4) 

Here, we consider the simplest n-dimensional spherically symmetric metric with one met- 
ric function, which may be written as 

ds^ = -f{v, r)dv^ + 2drdv + r^jdO'de^, (5) 
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where < r < oo is the radial coordinate, — oo < f < cxo is an advanced time coordinate, 
and 'jijdd^dd^ is the hne element of the (n — 2)-dimensional unit sphere. We want to intro- 
duce the general solutions of third order Lovelock gravity in the presence of a perfect fluid. 



In Refs. 



and 



111 ], the reverse has been done in 4-dimensional and ra- dimensional Ein- 



stein gravity, respectively. That is, a spherically symmetric solution with a metric function 
f{v, r) = 1 — m{v, r)/r is accepted and the properties of the energy-momentum tensor have 
been investigated. The energy-momentum tensor which we consider here is the sum of two 
noninteracting components namely the Vaidya null radiation and a type I perfect fluid given 
as m 



(6) 



where = (1, 0, 0, 0) and w^, = (//2, — 1, 0, 0) are two linearly independent future 



pointing null vectors in n dimensions which satisfy v^w^ 



-1. The stress-energy tensor 



has been chosen in such a way that T^yV^v 



and T^^w^'w" 



0", and therefore there 



is energy flux only along one of the null directions. It is of precisely the form which gives 
the charged Vaidya solution of Einstein gravity Jl3[| for p = p, and reduces to the energy- 
momentum tensor which gives the Vaidya metric [ij] for p = p = 0. The stress-energy tensor 
([6]) satisfies the dominant or weak energy conditions if the conditions [p > 0, —p < p < p 
and cr > 0] or [p > 0, p + p > and a > 0] are met, respectively. 

Using the field equation ([T]) in a unit system with a' = 1 and defining a'2 = 02/ {n — 
3)(n — 4) and a'^ = a^/S^n — 3)...(n — 6) for simplicity, the Q^'" and components reduce 
to: 



(^-2) 



1 + ^(1 



- 3) + -2 (1 - /) + (1 - /) 



klp{v, r) 



(7) 



in -2) 



2a 



"3 



2?^ I ^2 ^4 



(8) 



where the prime and the dot denote the derivatives with respect to the coordinates r and v, 
respectively. Also, one may note that the angular components of Eq. ([T]) may be written 

as 



p 



1 



d 



[n — 2)r"~'^ dr 



(r^-'p). 



(9) 



In order to solve the field equations, we define the energy function e{v,r) as 

op r 

e{v,r) = / piv^ry^'^^dr. (10) 

Now, substituting p{v, r) from Eq. ([7]) into Eq. (ITOl) and integrating the result with respect 
to r, one obtains: 



e(^;, r) = r"-^ ^ [1 - /(t;, r)] + a,r' [1 - /(t;, r)]^ + y [1 - /(^, 0]^ ^ . (11) 



Differentiating Eq. ffTTj) with respect to f and using Eq. ([8]), one finds: 



which shows that the functions p{v,r) and cr(f,r) are not arbitrary for our ansatz metric 

By solving Eq. ( ITTi) . the metric function is found to be: 

f{v,r) = |l - (^V^Pr^iv^ + k{v,r)^ ^ V 7^/^ (^7 + /i;^(^^,r') + k{v,r)^ , 

(13) 

where 



The solution introduced by Eqs. ([5]) and (|T4|) is a general spherically symmetric solution of 
third order Lovelock gravity with the ansatz metric ([5]) in the presence of a type II perfect 
fluid, where p(t>,r), a{v,r) and p{v,r) are related to each other according to Eqs. (jH]) 
and 0121) . This solution contains all the previous solutions of third order Lovelock gravity 

n Q 

introduced in the literature [4, |7| and contains some new exact solutions which will be 
discussed in the rest of the paper. 

III. STATIC SOLUTIONS FOR LINEAR EQUATION OF STATE: 

In this section, we find the static solutions of third order Lovelock gravity in the presence 
of a type I (a = 0) perfect fluid. Knowing the equation of state, and using Eq. ([9]), one may 
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obtain the density function p(r), and therefore the energy function e{r) exphcitly. For the 
hnear equation of state p = wp, Eq. ([9]) reduces to 



^ [r"-V(r)] +{n- 2)wr''-^pir) = 0, (15) 

with the solution 

where the integration constant A^ is positive in order to have the weak and dominant energy 
conditions. Using Eq. ffTOj) . the energy function e{r) may be obtained as 

A^ 

e(r) = m — - — ; ^ — — r^— -: w(n — 2) ^ 1, 

= m + A^lnr; w{n-2) = l. (17) 

For w = 1, with = A^/(n — 3), the function k{r) becomes 
The solution given by Eqs. (fT3l) and (fT8l) is the asymptotically flat static charged black hole 



of third order Lovelock gravity introduced in 

For w = —1 with the choice of \^ = {n — 1)A, the function k{r) reduces to 

fc r =-1 + — ^ + — I A + r . 19) 

^ ' 2al 2al V r^-V ^ ^ 

Due to the weak and dominant energy condition, A should be positive and therefore the 
solution given by Eqs. (fT3!) and (fT9|l presents an asymptotically dS uncharged solution of 
third order Lovelock gravity. 



A. Black hole for w = : 

The static case with w = gives a new asymptotically fiat solution of third order Lovelock 
gravity. The metric function /(r) is the solution of the following equation 

r"-^ [r^ [1 - /(r)] + a^r^ [1 - f{r)f + y [1 - f{r)f] = AV + m. (20) 

The solution of Eq. (l20l) is given in Eq. (fT3|) with the following A;(r): 

, , , Sets Sctq / m A^ \ , . 
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One may show that the above solution is asymptotically flat. In order to study the general 
structure of this solution, we flrst look for the curvature singularities. It is easy to show 
that the Kretschmann scalar R^uXkR'^'^^'^ diverges at r = 0, it is flnite for r 7^ and goes 
to zero as r — * 00. Thus, there is an essential singularity located at r = 0. Also, it is 
notable to mention that the Ricci scaler is flnite everywhere except at r = 0, and goes to 
zero as r ^ 00. The event horizon(s), if there exists any, is (are) located at the root(s) of 

grr = f ^r) = 0: 

r!^-^ + a2r+"^ + — r!^"^ - AV+ - m = 0. (22) 
3 

In order to flnd the number of positive real roots of Eq. ( 122|) . one should note that if the 
parameters m and A are chosen such that there exists a positive real root for /(rext) = 
f'{fext) = 0, then Eq. (!22|) has one real root. Differentiating Eq. (!20!) with respect to r and 
using /(rext) = /'('"ext) = 0, one obtains: 

(n - 3)r:-' + {n- 5)a2Cf + f - 7)Cf - = 0. (23) 

One may flnd a relation between mgxt and Agxt for the case that /(r) has only one real root 
by omitting Text between Eqs. fl23|) and fl22l) for r+ = Text- For n = 7, the relation between 
mext and Aext is found to be 

I / / \ 1/2 
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Aext = n ( ~6a2 + 6\/aj + 4^3 - 12mcxt 2^2 + \ + - 12mext i , 



which is real and positive provided m < a^/S. When m < 0:3/8, the solution given by Eqs. 
( IT3l) and ( l2Ti) presents a naked singularity if A < Aoxt, an extreme black hole for A = Aext, 
and a black hole with inner and outer horizons provided A > Aext- For m > 0:3/3, the 
solution is a black hole with one event horizon. To be more clear on this explanation, one 
may see the diagram of /(r) versus r for these four cases in Fig. [TJ 

The temperature may be obtained through the use of the deflnition of surface gravity. 
One obtains: 



y ^ />+) ^ in - 3)rr^ + 02(n - 5)rl-' + f{n - 7)rr^ - A^ 

47r 47r(r!^-=' + 2o2r!^"^ + r!^"^03) ' ^ ^ 



The entropy of asymptotically flat black holes of Lovelock gravity is [1 



O [("-l)/2] 
ZTT \ ^ 



S = — ^ pa', I d-~'x^~gC,^, (25) 
'^^ p=i 



FIG. 1: /(r) vs. r for n = 7, Q2 = 0.5, Q3 = 0.4, m = .1 < 03/8, A < Aext, A = Aext, A > Aext and 
m = .5 > 03/3 from up to down, respectively. 



where the integration is done on the {n — 2)-dimensional spacehke hypersurface of Kilhng 
horizon, g^^ is the induced metric on it, ^ is the determinant of g^i, and Ck is the kth order 
Lovelock Lagrangian of g^,^. It is a matter of calculation to obtain the entropy as 



+ \^ (n — 4)r^ (n — 6)r^ / 



(26) 



where Vn-2 is the volume of (n — 2)-dimensional unit sphere. Using Eqs. ([23]) and fl2^ . one 
may calculate the integral of J TdS 



n 



which is proportional to the geometrical mass m. 

Here, it is worth to consider a peculiar property of the above solutions which does not 
occur in the lower order Lovelock gravity. The solution with vanishing geometrical mass 
m = presents a black hole with two inner and outer horizons provided A > Aext, where 
Aext can be obtained through the use of Eqs. fl22|) and fl23l) as: 

Aelt = r:-,' {rl, + «2relt + f ) , (28) 



2 (^-6) J ^ , / 2 4(n-4)(n- 



The radius of the extreme black hole is zero for lower order Lovelock gravity and n = 8 as 
one may see from Eq. fl29p . and therefore the solution in these cases presents a black hole 
with only one horizon if A > 0. However, the solution for n = 7 with positive or n > 8 
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with negative as presents a black hole with two inner and outer horizons if A > Acxt- This 
is a peculiar feature of third order Lovelock gravity which does not occur in the lower order 
Lovelock gravity in the presence of a perfect fluid. 



B. Black hole for w = {n — 2) ^ : 

As in the case of w = 0, the static case with w = {n — 2)~^ gives a new solution of third 
order Lovelock gravity with 

, 1 3 I/O 3ao(m + A^lnr) 
k(v, r) = - + -Y' H — 5 (30) 

which is asymptotically flat. Again, the Kretschmann scalar R^^xi^R^^^'^ diverges at r = 0, 
it is finite for r 7^ and goes to zero as r ^ 00. Thus, there is an essential singularity 
located at r = 0. The event horizon(s), if there exists any, is (are) located at the root(s) of: 

r+"^ + a2r+"^ + — r'^"^ - m - A^ In r+ = 0. (31) 
3 

Following our leading procedure in the above subsection for n = 7, the relation between 
mext and Aext is found to be 

The solution given by Eqs. f|T3l) and fl30l) presents a naked singularity if m < rriext, an 
extreme black hole for m = mext, and a black hole with inner and outer horizons provided 
m > rriexf The diagrams of /(r) versus r for these cases are shown in Fig. [21 Again, one 
may note that J TdS is proportional to the geometrical mass of the black hole. 



C. Strange quark matter: 

For completeness, we consider the linear equation of state of a strange quark matter 
(SQM). This is due to the fact that near the singularity, matter may be in the highest known 
density form, which brings in the SQM. In 4 dimensions, the SQM fluid is characterized by 
the equation of state p = w{p — 45/3) where B is the bag constant indicating the difference 



between the energy density of the perturbative and nonperturbative QCD vacuum 16|. 



Here, we consider the equation of state p = w{p — 4i?) as a generalization of equation of 
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FIG. 2: /(r) vs. r for n = 8, 02 = 0.5, Q3 = 0.3, A = .5, m < rriext, ?7i = rrtext, and m > mext from 
up to down, respectively. 



state of SQM in n dimensions, where 5 is a constant [l7|, [l8|. In this case, the functions 

p(r) and e{r-) reduce to: 



e{r) = m + 
= m + 



p{r) 



+ 



B 



(32) 



{l + w){n-l) [(n - 2)u; - l]r("-2) 
+ A'' Inr; 



(1 + w){n - 1) 

respectively. The function k{r) for the above two cases becomes 



w{n- 2)^1, 
{n - 2)w = 1, 



(33) 



k{r) 



3^3 



QwBal 



3ao / m , , , , 

+ 7r4{—^+ , ...u^ 1 ; w{n- 2)^1, 



2ci\ a\{\^w){n-\) V^""^ ^(i+«')(n-2) 

3a| (m -I- A^ Inr) 



3^3 6w-Bao 



+ 



2a|r"-i 



n — 2^w = 1. 



(34) 



2a| a|(l-Fw)(n-l) 

The i?-term in the above expressions guarantees that the solution is asymptotically (A)dS 
even in the cases of w = 1, w = and w = {n — 2)^^, which are asymptotically flat for 
normal matter. 



IV. RADIATING SOLUTIONS: 

Now, we want to find the radiating solutions of third order Lovelock gravity in the 
presence of energy-momentum tensor ([6]) for linear equation of state with a{v, r) given in 



10 



Eq. ([8]). In this case, the energy density e{v,r) becomes 

eiv, r) = miv) — — ; — — — -— -: win — 2) 7^ 1. 

= m(f) + A^(f)lnr; ^(n - 2) = 1. (35) 

When w = — 1, the solution is the asymptotically (A)dS uncharged Vaidya-type solution 
introduced in Ref. 

For w = 1 with q^{v) = \^{v)/{n — 3), the solution is a new solution which may be 
called the charged Vaidya-type solution of third order Lovelock gravity for which the metric 
function is given by Eq. (fT3l) with: 

In this case, the energy-momentum tensor is the sum of av^Vy and the Maxwell energy- 
momentum tensor 

T^r) = F,yF^ ' - \g,.F,^FP\ (36) 

for a point charge g(f ) with electromagnetic field tensor F^i, = d^A^ — dyAfj, and the following 
potential: 

A.. = ^Hs',. 

Using Eq. f|T2l) . one obtains 



2/t2 \.fn-2 ^2n-5 ' ' ^ ' 



which shows that when q < 0, the weak and dominated energy conditions will be satisfied. 
While for q >0, the energy conditions is satisfied for r > r^, where = (2gg/m)^/*^""'^\ But 
in realistic situations the particle cannot get into the region r < because of the Lorentz 
force, and therefore the weak and dominant energy conditions are satisfied for the charged 



Vaidya-type solution [191]. In the limit of 02 = as = 0, this solution reduces to the charged 
Vaidya solution of Einstein gravity introduced in Ref. [3]. 

The nature of the singularity (to be naked or hidden) can be characterized by the existence 
of radial null geodesies coming out of the singularity. The nature of the singularity is exactly 
the same as the uncharged solutions explained in Ref. 7(]. Thus, we only compare the 
strength of the singularity with the case of uncharged solution. The solution satisfies a 
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strong curvature condition (SCC) [20] or limiting focusing condition (LFC) [2l|] provided 
the limits of or r$ are positive, respectively, where r is an affine parameter and $ is 



Using the fact that dr /dr = {dv/dT)f/2, one can show that 



2{n - 2)f (dr 



2 



and the radial null geodesic satisfies the differential equation 

d^r 2/ f dr^^ 



(39) 



dX^ P \dX, 

Now, we consider the strength of the singularity for the following two cases: 

1. m{v) = mo^(f)f"~^ and q^{v) = q^O {v)v'^^'^~^\ where mo and go are two arbitrary 
constants, and 6{v) is the step function which is 1 for f > and zero for f < 0: In this case, 
the limits of / and / are 1 and zero as r = f — > 0, respectively. Thus, one finds that the 
limit of r$ and r^<l> are zero as r goes to zero, and therefore neither the SCC nor LFC are 
.a«*d ^on. a .ad.al geode.c. The unchanged solution .a«.fle= .he LFC Q. wMle 
the charged solution does not satisfy LFC, and therefore charge weakens the strength of the 
singularity. 

2. m(f ) = mo^^(f)f"~^ and q'^{v) = go^(f )f where again mo and go are two arbitrary 
constants: In this case, the limits of / and / as r and v go to zero are: 

lim/ = 1 

W ^ _ (3aI)V3|(„-4)^„-(2„ 7),g]) ^ ^ 

r^o 3 (mo — go) ' 

Using Eqs. ( l39l) and ( HOl) . one may show that the radial null geodesic near r = t> = is 
r ~ (/o)~^ln(r + 1). Using this result and Eqs. ( |38l) and ( HOl) . one finds that the limit of 
r$ is positive while the limit of is zero as r goes to zero, and therefore as in the case 
of uncharged solution [3], only LFC is satisfied along a radial null geodesic. Although the 
radiating solution with w = is a new solution with 

^ ^ 2al 2a\ \ r"-i r"-2 

but similar calculations show that the strength of the singularity for this case is exactly the 
same as when \{y) = ^], and therefore we will not present them here. 
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As in the case of static solutions, the effect of using the equation of state of SQM instead 
of normal matter is to make the solutions asymptotically (A)dS. These solutions reduce to 
the radiating solutions of Einstein gravity in the presence of SQM introduced in Ref. 18 1. 



V. CLOSING REMARKS 

We considered the third order Lovelock gravity for a spherically symmetric spacetime in 
the presence of a type II perfect fluid. Due to the fact that the ansatz metric ([5]) had only 
one unknown function, the nonvanishing components of Qfj,^ of the LHS of Eq. ([1]) were 
related to each other according to the following equations: 

1 (9 

(n-2)r"-3 5r^ ^ ' 

Thus in our analysis, the functions p, p and a in the energy-momentum tensor ([6]) were 
not arbitrary. Nevertheless, we found a general solution for linear equation of state which 
contains all the known solutions of third order Lovelock gravity, and also some new static 
and radiating solutions. We investigated the properties of static solutions for w = and 
w = {n — 2)"^ with zero and nonzero B, and found that these solutions may be interpreted 
as black holes with inner and outer horizons, extreme black holes, naked singularities or 
black holes with one horizon provided the metric parameters are chosen suitable. We found 
that the solution with w = and m = may present a black hole with inner and outer 
horizons. This is a peculiar feature of third order Lovelock gravity which does not occur 
in lower order Lovelock gravity in the presence of perfect fluid. We also presented the new 
radiating solutions for w = and w = (n — 2)^^, and compared the singularity strengths 
of these solutions and the known radiating solutions of third order Lovelock gravity in the 
literature. Although naked singularity is inevitably formed in third order Lovelock gravity 
for the general energy- momentum tensor ([6]), the strength of the singularity depends on the 
rate of increase of e{v, r) with respect to v. For the case of charged solutions when the rate 
of increase of charge is large enough, then the strength of the singularity is weaker than 
that of the uncharged case. That is, the charge with a suitable rate of increase weakens 
the strength of the singularity. We also found that the presence of SQM instead of normal 
matter changes the asymptotic behavior of the solutions from flat to (A)dS. 
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Although we investigated the solutions ior w — 0, w — 1 and w — (n — 2)~^, one 
may consider the solutions for other values oi w. If one likes to obtain the most general 
spherically symmetric solution of Lovelock gravity in the presence of a type II perfect fluid 
with arbitrary p, p and a, then one should consider the following metric: 

ds'^ = -A{v, r)dv^ + 2B{v, r)drdv + r^C(v, r)-fijd9'9^ . 
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